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ABSTRACT 



Context. The ionization fraction (i.e., the electron abundance) plays a key role in the chemistry and dynamics of molecular clouds. 
I ■ Aims. We study the H'^'CO^, DCO^ and HOC^ line emission towards the Horsehead, from the shielded core to the UV irradiated 

O ' cloud edge, i.e., the Photodissociation Region (PDR), as a template to investigate the ionization fraction gradient in molecular clouds. 

Methods. We analyze an IRAM Plateau de Bure Interferometer map of the H'^^CO^ 7=1-0 line at a 6.8" x 4.7" resolution, comple- 
C/3 \ mented with mAM-30m H''CO+ and DCO+ higher-y line maps and new HOC^ and CO^ observations. We compare self-consistently 

. the observed spatial distribution and line intensities with detailed depth-dependent predictions of a PDR model coupled with a non- 

local radiative transfer calculation. The chemical network includes deuterated species, '^C fractionation reactions and HCO^/HOC* 
isomerization reactions. The role of neutral and charged PAHs in the cloud chemistry and ionization balance is investigated. 
Results. The detection of HOC^ reactive ion towards the Horsehead PDR proves the high ionization fraction of the outer UV irra- 
diated regions, where we derive a low [HCO^]/[HOC^]-75-200 abundance ratio. In the absence of PAHs, we reproduce the obser- 
vations with gas-phase metal abundances, [Fe+Mg+...], lower than 4x10"' (with respect to H) and a cosmic-rays ionization rate of 
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Conclusions. The ionization fraction in the Horsehead edge follows a steep gradient, with a scale length of ~0.05 pc (or ~25"), from 
[e^l-lO""* (or ne~l-5 cm"^) in the PDR to a few times ~10"' in the core. PAH" anions play a role in the charge balance of the cold 
^Nj ' and neutral gas if substantial amounts of free PAHs are present ([PAH]>10"*'). 
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1. Introduction Rotational line emission of molecular ions such as DCO^ 

and HCO^ have been traditionally used to estimate the 
^ , The electron abundance ([e ]=n,/nH) plays a fundamental ionization fraction in molecular clouds because (i) they 
d ■ role in the chemistry and dynamics of interstellar gas. The abundant and easily observable (//) dissociative recom- 

degree of ionization determines the preponderance of ion- bination is their main destruction route, and thus their 
neutral reactions, i.e., the main form ation route for most chem- abundances are roughly inversely proportional to the elec 
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ical species in molecular clou ds ( ,Herbst & Klemperer 1973| tron abundance (e.g., LGuelin et al. 1982) IW ootten et al. 19821 
Oppenheimer& Dalgarno 1974^ . In addition, the ionization |de Boisanger et al. 1996 Williams et al 1998, Caselh et al. 



fraction constrains the coupling of matter and magnetic fields, 1998, | Maret & Bergin 2007] iHezareh et al. 2008 ). On the other 



which drives the dissipation of turbulence and the transfer of an- ^and, the presence of reactive ions (species such as HOC+ or 

gular momentum, thus having crucial imp hcations in protostellar ^0+ that react rapidly with H2) is predicted to be a sensitive in- 

collapse and accretion disks (e.g., .Balbus & Hawley 1991j . ^icator of high ioniz ation fraction regions, e.g., theUV irradiated 

High-angular resolution observations of interstellar clouds cloud surfaces (e.g- ISmith et~2 002 Fue nte et al. 20 03 ). 

reveal steep density, temperature and turbulence gradients as ^ , ..... ^ . ,■ n 

well as sharp chemical variations. Accordingly, the electron ]" "I'^J^i" constrain the ionization fraction gradient from 

abundance should vary within a cloud depending on the relative "^"^^f *^ '^^^^'^ chemistry and physics can not be simplified 

ionizing sources and prevailing chemistry. ^ the charge balance depends on parameters such as 

the penetration of UV radiation, the cosmic -rays ionization rate 

(O and the abundance of key species (e.g., metals and PAH). 

Send offprint requests to: e-mail: goicoechea@damir . iem. csic. es 

* Based on observations obtained with the IRAM Plateau de Bure Compared to other works, in this paper we determine the 

interferometer and 30 m telescope. IRAM is supported by INSU/CNRS ionization fraction gradient by direct comparison of H'^CO+ 

(France), MPG (Germany), and IGN (Spain). and DCO^ high-angular resolution maps and HOC^ pointed ob- 
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Table 1. Observation parameters of the PdBI maps shown in Figure[T] 



Molecule 


Transition 


Frequency 


InstiTjment 


Config. 


Beam 


PA 


Vel. Resol. 


Int. Time" 


T 

^ sys 


Noise*'^ 


Obs. date 






GHz 






arcsec 




km s"' 


hours 


K 


K 






1-0 


86.754288 


PdBI 


C&D 


6.8x4.7 


13 


0.2 


6.5 


150 


0.10 


2006-07 


HCO 


lo.i 3/2,2 -Oo.o 1/2,1 


86.670760 


PdBI 


C&D 


6.7 X 4.4 


16 


0.2 


6.5 


150 


0.09 


2006-07 



" We observed a 7-field mosaic centered on the IR peak at ajooo = 05''40'"54.27% ^2000 = -02°28'00" (Abergel et al. 2003) with the following 
offsets: (-5.5", -22.0"), (5.5", -22.0"), (11.0", 0.0"), (0.0", 0.0"), (-11.0", 0.0"), (-5.5", 22.0") and (5.5", 22.0"). The total field-of-view is 
80.1" X 102.1" and the half power primary beam is 58.1". The mosaic was Nyquist sampled in declination at 3.4 mm and largely oversampled in 
right ascension. This maximizes the field of view along the PDR edge while the oversampling in the perpendicular direction eases the deconvolu- 
tion. On-source time computed as if the source were always observed with 6 antennae. * The noise values refer to the mosaic phase center (mosaic 
noise is inhomogeneous due to primary beam correction; it steeply increases at the mosaic edges). 

Table 2. Observation parameters of the IRAM-30m observations. 



Molecule 


Transition 


Frequency 


Instrument 




Beff 


Resol. 


Resol. 


Int. Time 


Noise ' 


Observing 


Obs. date 






GHz 








arcsec 


km s"' 


hours 


K 


Mode 




HCO+ 


7=1-0 


89.188523 


30m/A100 


0.95 


0.78 


27.6" 


0.20 


4.7 


0.02 


ON-OFF 


2008 


HOC+ 


7=1-0 


89.487414 


30m/A100 


0.95 


0.78 


27.5" 


0.20 


4.7 


0.02 


ON-OFF 


2008 


CO+ 


2, 5/2-1,3/2 


236.062578 


30m/A230 


0.91 


0.52 


10.4" 


0.20 


4.7 


0.05 


ON-OFF 


2008 


H'3C0+ 


7=1-0 


86.754288 


30m/AB100 


0.95 


0.78 


28.4" 


0.20 


2.6 


0.10 


OTF map 


2006-07 


h13C0+ 


7=3-2 


260.255339 


30m/HERA 


0.90 


0.46 


13.5" 


0.20 


5.9 


0.06 


OTF map 


2006 


DCO+ 


7=2-1 


144.077289 


30m/CD150 


0.93 


0.69 


18.0" 


0.08 


5.9 


0.18 


OTF map 


2006 


DCO+ 


7=3-2 


216.112582 


30m/HERA 


0.90 


0.52 


11.4" 


0.11 


1.5 


0.10 


OTF map 


2006 



The noise (in T^b scale) refers to the channel spacing obtained by averaging adjacent channels to the velocity resolution given in the tables. 



servations, with detailed depth-dependent chemical and radia- 
tive transfer models covering a broad range of cloud physical 
conditions. Indeed, the observed field-of-view contains the fa- 
mous Horsehead PDR (the UV illuminated edge of the cloud) 
and a dense and cold core discovered by us from its intense 
DCO^ line emission (Pety et al. 2007'). Due to its simple ge- 
ometry and moderate distance (c/^400pc), the Horsehead PDR 
and associated core are good templates to study the steep gradi- 



ents expected in molecular clouds (e.g., Pety et al. 2005| 2007, 



IGoicoechea et al. 2006IIGerin et al. 2009I I. 



The paper is organized as follows. The observations are pre- 
sented in Sect|2] and the models used to interpret them are de- 
scribed in Sect. [3] The chemistry of H'^CO+, DCO+ and HOC+ 
(our observational probes of the ionization fraction) is analyzed 
in Sect.|4] In Sect. |5] we investigate the role of metals, PAHs and 
( on the electron abundance determination. The main results and 
constrains are presented in Sect.|6]and discussed in Sect|7] 

2. Observations 

2.1. Observations and data reduction 

Tables [T]and|2]summarize the observation parameters of the data 
obtained with the PdBI and the IRAM-30m telescope that we 
shall study in this work. The H'^CO^ 7=1-0 line emission map 
was first presented in Gerin et al. (2009). Frequency-switched, 
on-the-fly maps (OTF) obtained at the IRAM-30m were used to 
produce the short-spacings needed to complement a 7-field mo- 
saic acquired with the 6 PdBI antennae in the CD configuration 
(baseline lengths from 24 to 176 m). Correlator backends were 
used (VESPA for IRAM-30m observations). The high angular 
resolution PdBI H'^CO^7=l-0 map complements our previous 
H'^CO+y=3-2 and DCO+ 7=2-1 and 3-2 maps taken with the 
IRAM-30m telescope and first presented in Pety et al. (2007). 

In this work we present new IRAM-30m deeper integrations 
in the HOC+, H'^CO+and HCO+ 7=1-0 lines, and an upper limit 
for the CO^ emission towards the PDR (defined here as the HCO 
emission peak; IGerin et al. 20091 1 . The position switching ob- 
serving mode was used. The on-off" cycle duration was 1 minute 



and the off'-position offsets were (Aa, A6) - (-100", 0"), i.e., the 
Hii region ionized by crOri and free of molecular gas emis- 
sion. Position accuracy is estimated to be ~3" for the 30m 
data and better than 0.5" for the PdBI data. The data process- 
ing was done with the GILDA^ softwares (e.g., |Pety 2005b) l. 
The IRAM-30m data were first calibrated to the T^ scale using 
the chopper wheel method (Penzias & Burr us I973I I. and finally 
converted to main beam temperatures T^b using the forward 
and main beam efficiencies Feft and Z?eft displayed in Table |2] 
(e.g.. lGreve et al. I9981 l. The amplitude accuracy for heterodyne 
observations with the IRAM-30m telescope is ~IQ%. PdBI data 
and short-spacing data were merged before imaging and decon- 
volution of the mosaic, using standard techniques of GILDAS and 
used in our previous works (see e.g., |Pety et al. 2005) . 



2.2. DCO^ andH^^CO^ spatiai distribution, HOC^ detection 

Figure m shows Hi^CO^ 7=1-0, HCO lo.i-Oo.o (PdBI) and 
DCO ^ 7=2-1, 3-2 i ntegrated line intensity maps (IRAM- 
30m; Pety et al. 2007 1, as well as the Aromatic Infrared Band 
emission (AIB, observed with ISOCAM, [Abergel et al. 2003| l 
that traces the UV illuminated edge of the cloud, i.e., the 
PDR. The DCO^ emission is concentrated in a narrow, 
arclike structure of dense and cold gas behind the PDR 



(Petyetal. 2007 1. Hence, it shows a very different spatial dis- 



tribution than the emission of "PDR tracers" such as C2H, C4H, 
C-C3H2 d Pety et al. 2005|), HCO r adicals ( Gerin et al . 2009'), vi- 
brationally excited H? (IHabart et al. 2005) or the AIB emission 
(Compiegne et al. 2008). The Hi^^CO+ 7=1-0 emission fol- 
lows the DCO^ distribution and it mostly delineates the dense 
core that coincides with the DCO^ emission peak. Nevertheless, 
while DCO^ is not detected in the illuminated edge, H'^CO^ 
does show a faint emission in the PDR. Therefore, the small 
field-of-view shown in Fig. [1] contains two different environ- 
ments: a warm PDR and a cold core shielded from the external 
UV radiation field. In the following sections we analyze these 



See http : //www . iram. fr/IRAMFR/GILDAS 
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Fig.l. DCO+ 7=3-2 and 2-1 (IRAM-30m; from |Pety et al. 2007| ) H'^CO^ 7=1-0 (P dBI) and 3-2 (IRAM-3 0m) line integrated 
intensity maps. Aromatic Infrared Band emission (ISOCAM, from Abergel et al. 2003 ) and HCO (PdBI, from lGerin et al. 20091 1. 
Maps have been rotated by 14° counter-clockwise around the projection center, located at {5x, 6y) = (20", 0") to bring the illumi- 
nated star direction in the horizontal direction. The horizontal zero has been set at the cloud edge {6x-0"). The H'^CO^, DCO^ and 
HCO emission is integrated between 10.1 and 11.1 kms Integrated intensities are expressed in the T,„i scale. Contour levels are 
displayed on the grey scale lookup tables. The red vertical line shows the PDR edge and the green crosses shows two representative 
positions: the "shielded core" (the DCO^ emission peak at ^x~45"; |Pety et al. 2007 1 and the "PDR" (the HCO emission peak at 
(5x~15"; ,Gerin et al. 2009J . The dashed blue line shows the horizontal cut analyzed in this work. 



emission maps to determine the ionization fraction gradient in 
the region. 

Figure |2] shows long integration spectra of the HOC^, 
H'^CO+ and HCO+ 7=1-0 lines towards the PDR. This 
is the first detection of HOC^ reactive ion towards the 
Horsehead, and adds to previous detections in interstellar en- 
vironments with high electron abundances d Woods et al. 19831 
Ziurys & Apponi 1995[ [Fuente et al. 20 0l IRizzo et al. 2"003l 
Savage & Ziurys 2004| ILiszt et al. 20041). H^^CO^ lines are op- 



Table 3. Main spectroscopic parameters of the studied lines. 



tically thick, as shown by the low U^^CO^fH^^CO^ 7=1-0 line 
intensity ratio (~7), much lower than the expected ^^C/^^C^60 
abundance ratio (Langer & Penzias 1990 Savage et al. 2002| 
and references therein). The large opacity of H^'^CO^ lines even 
towards the PDR justifies the use of H'^CO^ lines as tracers of 
the HCO^ abundance. 



Species 


Transition 


Frequency 


A.J 




^upp 




7upp ~ 7]ow 


(GHz) 


(s-') 




(K) 


HCO+ 


1-0 


89.188523 


4.2x10" 


-5 


4.3 


HOC+ 


1-0 


89.487414 


2.2x10- 


-5 


4.3 


H'''CO+ 


1-0 


86.754288 


3.9x10- 


-5 


4.2 


C0+ 


2(5/2)-l(3/2) 


236.062578 


4.7x10" 


-4 


17.2 


H'''CO+ 


3-2 


260.255339 


1.3x10" 


-3 


25.0 


DCO+ 


2-1 


144.077289 


2.1x10" 


-4 


10.4 


DCO+ 


3-2 


216.112582 


7.7x10" 


-4 


20.7 



3. Analysis: Models 

In this work we couple the depth-dependent abundances pre- 
dicted by a PDR model (for the varying physical conditions pre- 
vailing in the Horsehead edge) with detailed excitation and ra- 
diative transfer calculations adapted to the cloud geometry. This 
technique allows us to analyze different chemical models by di- 
rect comparison with observed line intensities. This methodol- 
ogy was introduced to study our interferometric CS and C'^O 



maps of the Horsehead edge jGoicoechea et al. 20061) . It enables 
to observationally benchmark the abundance gradients predicted 
by chemical models, even if it does not produce perfect fits to 
line profiles in all cloud positions. In this paper, we analyze a 
horizontal cut of the H'^CO^ and DCO^ line emission along the 
direction of the illuminating star (6y=l5"). Figure [1] shows that 
this cut (blue dashed line) goes across the DCO^ emission peak 
{6x ~45"), which we identify as the "shielded core", and across 
the HCO emission peak, the "PDR" (6x~l5"). 
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Fig.2. HOC+ and H'^CO+ 7=1-0 lines towards the Horsehead 
PDR (upper and middle panels) observed with the IRAM- 
30m telescope. Solid lines are radiative transfer mod- 
els with Ti=60K, n(H2)=5xlO'*cm ^ n(H)=500cm-^ and 
[e ]=5xlO"^. Three different abundances are shown, grey thick 
line: [HOC+]=4.0xlO-i2 and [H'^CO+l^l.SxlO "; red dashed 
line: abundances x2; blue thin line: abundances -h2. For com- 
pleteness, the HCO^ 7=1-0 line towards the PDR is also shown 
(lower panel). This transition is very opaque, as shown by the 
low H'2C0VH'^C0+ 7=1-0 hne intensity ratio (~7). The result- 
ing line profile is thus broadened and it suffers from scattering 
by low-density foreground gas that we do not model here. 



3.1. Geometry and density gradient 

The Horsehead edge has an almost edge-on geometry with a 
line-of-sight depth of Idepth -0. 1 pc (e.g., Habart et al. 2005 ) and 
a spatial scale in the plane of the sky of ^0.002 pcarcsec 
We determine the density profile from observations by fitting 
the 1.2 mm dust continuum emission (IRAM-30m/MAMBO) 
along the 6y-l5" direction ( jHiTy-Bla nt et al. 2005). In this fit, 
we adopt a dust opacity per unit (gasH-dust) mass column density 
of /fi.2=0.003 cm^ g"' at 1.2 mm (computed for "MRN grains": 
Mathis, Rumpl & Nordsieck 1977, see below), our best knowl- 
edge of the dust grains temperature (from ~15K in the core 
to -30 K in the PDR; 



e.g., Ward-Thompson et al. 2006| ) and a 



power-law density profile nH(r)= n(H)+2n(H2) r~^, where r 
is the distance from the shielded core towards the illuminated 
edge of the cloud. Best fits are obtained for a steep density gra- 
dient in the cloud edge (p ^3) and a flatter one towards the core 
(p ^0.5). The turnover point occurs at core radius of r ^0.04 pc 
(or 6x ^23" in the maps). The resulting density gradient used 
in the photochemical and radiative transfer models is shown in 
Fig. m In the next sections we constrain the ionization frac- 
tion gradient in the cloud by comparing synthetic and observed 
H'^CO^ and DCO^ spectra along the same cut. 



3.2. Piiotocliemical models 

We have updated the Meudon PDR code to model our obser- 
vations of the Horsehead. The code has been described in de- 
teil^elsewhere (e.g., Le Bourlot etal. 19931 |Le Petit et al. 20061 
[Goicoechea & Le Bourlot 2007) and benchmarked against other 
PDR codes by Rollig et al. (2007). In this section we summarize 
the most relevant upgrades and model features for this work. 



3.2.1 . UV radiative transfer and dust properties 

The code solves the UV radiative transfer taking into account 
dust scattering and gas absorption. Anisotropic scattering of UV 
photons by dust grains is included by explicity calculating the 
wavelength-dependent grain albedo and ^-asymmetry param- 
eters jGoicoechea & Le Bourlot 2007i This enables the spe- 
cific computation of the UV radiation field (continuum+lines) 
and thus, the direct integration of consistent photoionization 
and photodissociation rates. We use two types of dust popu- 
lations: (0 a mixture of graphite-Hsilicate grains and (//) PAHs 
(see next paragraph). More precisely, we adopt a power-law size 
distribution (n(a) oc a^^'^) with minimum and maximum radius 
of ~5 and ~250nm respectively (for graphite-i-silicate grains). 
Wavelength-dependent optical properties (Q efficiencies and g 
factors) are interpolated from Laor & Draine (1993) tabulations. 
With a standard gas-to-dust mass ratio (~100), this grain mixture 
("MRN grains") reproduces the main characteristics of the stan- 
dard interstellar extinction curve with A^hMv =1-9x10^' cm"^ 
and Rv =3.1. 

In order to complete our description of the dust popula- 
tions, in this work we have also added smaller aromatic grains. 
Observationally, the AIB emission towards the Horsehead (pro- 
duced by free PAHs according to the most accepted the- 
ory; |LSgeF&^geFT984 lAUamandola et al. 19851 ) clearly sep- 
arates the Hii region and PDR (where the emission is bright) 
from the regions shielded from UV radiation, where no AIB 
emission is detected ( Abergel et al. 2003 IHabart et al. 20051 
Compiegne et al. 2007 120081 ). However, the size distribution 
and PAH abundance in dense regions shielded from UV ra- 
diation are uncertain. It may vary from "negligible", if PAHs 
coagulate into larger PAH aggregates ( [Boulanger et al. 1990 



Rapacioli et al. 2006) to "high" abundances (though they wil 
not be detected in the mid-IR due to the lack of UV photons to 
excite them). We used the following PAH properties: a size dis- 
tribution with ~0.4 and ~1.2nm radii limits (Desert et al. 1990] ) 
and optical properties from Li & Draine (2001). This size 
distribution is compatible with PAHs having a mean radius 
of ~0.6nm and A^c~100 carbon atoms assuming A^c-500a^ 
dBakes & Tielens 19941 ). The extinction curve and the efficiency 
of the photoelectric heating mechanism depend on the mass frac- 
tion put into PAHs dBakes & Tielens 19941 ). Depending on the 
PAH abundance, they contribute to the total dust mass by ~1% 
for [PAH] = 10-^ and -10% for [PAH] = 10-''. 

3.2.2. Cliemical network and elemental abundances 

Once the UV field is determined at every cloud position, 
steady-state chemical abundances are computed for a given 
network of chemical reactions. The model also computes 
the temperature profile by solving the thermal balance be- 
tween the most important gas heating and cooling mechanisms 
dLe Petit et al. 2006 ). Our chemical network contains ~160 
species and ~2000 reactions. It includes deuteration, '^C frac- 
tionation dOraedel et al. 1982| i and HCO+/HOC+ isomerization 
reactions. When available, we used the photodissociation rates 
given by van Dishoeck (1988), which are explicitly calculated 
for the Draine's interstellar radiation field (ISRF). The most crit- 
ical reaction rates for our determination of the ionization fraction 
are listed in Table |4] Most reactions were checked against OSU 
(E. Herbst and co-workers) and UDFA dWoodall et al. 20061 net- 
works. Besides, we benchmarked our network with more ex- 
tended ones by comparing the predicted abundances of simple 
species such as CO and DCO^. 
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Following Flower & Pineau des Forets (2003), we have 
also included interactions (~60 reactions) of gas phase species 
with very small aromatic grains (neutral PAH and singly 
charged PAH'^). In particular, we take into account PAH- 
gas processes such as neutralization reactions of atomic and 
molecular cations on PAH , PAH electron attachment and 
photodetachment of PAH" and PAH by UV photons. Such 
processes can play a significant role in the ionization bal- 



Table 5. Standard conditions and gas-phase elemental abun- 
dances. Abundances, noted [x], refer to H. 



ance of dense molecu l ar clouds (e.g., Lepp & Dalgarno 1988 
IBakes & Tielens 19981 iFloweretal. 20071 Wakelam & Herbst 
2008, IWolfire et al. 2008]) . We have not included larger grains 
in the network in order to isolate the role of PAHs in the gas- 
phase chemistry. We thus assume that recombinations of ions 
with grains are much less frequent than recombinations with 
electrons and PAH . This is partially justified by the fact that, 
according to their size and mass, the fractional abundance of 
"MRN grains" is low: ng/ne ~10""^ and their effective cross sec- 
tion per H nucleus is (ng/«H)7ra^ ~10"^' cm"^ (the product of 
the grains abundance and the mean grain cross section). Grain 
growth towards the denser parts of the cloud will result in even 
lower grain abundances and smaller effective cross sections if the 
gas-to-dust mass ratio has to be preserved: e.g., (wg/nn) ~10"'^ 
and (ng/nii)7Ta^ xilO"^^ cm"^ if the grain radii a increase by ~10. 
Therefore, the resulting lower abundance of charged grains and 
their smaller effective cross section for ion-grain recombinations 
will not alter our estimations of the ionization fraction much. 



Table 4. Key chemical reaction rates adopted in this work. 



Reaction 



Rate [cm^^ s"'] 



HCO+ + e- ^ CO + H 2.4 X lO-'OOOK/T)"^"" 

HCO+ + PAH- ^ PAH + CO + H 1.4 x 10-'*(300 K/T)" * 

M+ + e- ^ M -H hv 3.7 X 10-'-(300 K/T)"*' 

M+ + PAH- ->M + PAH 1.0 X 10-'*(300 K/T)" * 

C+ + H2O -> HCO+ + H 8.9 X 10-'°(300 K/T)" 

C+ + H2O ^ HOC+ + H 1 .8 X 10-''(300 K/T)"^" 

CO+ + H2 ^ HCO+ + H 7.5 X 10-'° 

CO+ + H2 ^ HOC+ -F H 7.5 X 10-'° 

HOC+ + H2 -» HCO+ + Ha 3.8 x 10-'°" 

' Rates are from the UDFA {U MIST Database for Astrochemistry) and 

OSU (Ohio State University) data bases unless indicated. 

" Cited in the te xt as a(HCO'^), it also applies to DCO+ and H'3C0+. 

* Rate is from Flower & Pinea u des Forets 20031 

' Rate is from .Smith et al. 2002] 



The adopted elemental abundances are shown in Table |5| 
Low ionization potential heavy metals (<8 eV; Fe^, Mg^ or Na^) 
are all represented by a single element, "M^". In our model, 
such metals slowly recombine with electrons (through radiative 
recombinations), can be photoionized and can exchange charge 
with ions and neutrals (including PAHs). However, they are as- 
sumed to be chemically inert and thus do not form metallic 
molecules (e.g., |Qppenheimer & Dalgarno 1974[ |. 

Once the physical and geometrical parameters of the cloud 
are constrained, the only free parameters in the model are the 
cosmic-rays ionization rate and the metal and PAH abundances. 

3.3. From abundances to spectra: mm radiative transfer 

We use the PDR model predictions (molecular abundance, 
n{H2), n(H), gas temperature and ionization fraction gradients) 
as input for a nonlocal radiative transfer calculation able to com- 



Parameter 



Value 



Radiation field 

Density nH(r)=n(li) + 2«(H2) 

Line of sight depth l^epih 

[He]=n(He)/nH 

[O] 

['^C] 

[N] 

[D] 

[S] 

['3C]=['2C]/60 
[PAH] 

[metals] = [M] = [Fe + Mg + ...] 
Cosmic ray ionization rate ( 



60 (Draine units) 
oc r--\ up to ~ 2 X 10'' cm-' 
0.1 pc 
1.00 X 10-' 
3.02 X 10-* 
1.38 X 10-" 
7.95 X 10-5 
1.60 X 10-"^ 
3.50 X 10-* 
2.30 X lO-*- 
variable: 0-10-^ 
variable: 10-" -10-^ 



pute DCO^ and H'^CO^ line intensities as a function of cloud 
position. Our radiative transfer code handles edge-on plane- 
parallel geometry, and accounts for line trapping, coUisional 
excitatioifl and radiative excitation by absorption of contin- 
uum photons. After the level populations are determined in each 
modeled slab, emergent line intensities along each line of sight 
are computed and convolved with the telescope angular resolu- 
tion at each frequency. A more detailed description is given in 
Goicoechea et al. (2006; Appendix). Since typical densities in 
the Horsehead (^lO'^-lO^ cm-^) are below the critical densities 
of the observed high-dipole moment molecular ions (a few ~ 10^- 
10'' cm-^ for the studied transitions) our approach allows us to 
properly take into account non-LTE excitation effects (e.g., sub- 
thermal excitation), as well as opacity and line profile formation. 

4. Chemistry of the ionization fraction probes 

4.1. H'^CO+ and DCO+ chemistry in tiie UV shielded core 

The detection of very bright DCO^ emission towards the 
shielded core ( Pety et al. 2007) implies cold gas temperatures 
(T;t ^10-20 K) and thus efficient HCO^ deuterium fractiona- 
tion (i.e., [DCO^]/[HCO+] »D/H). From our observations we 
infer [DCO+] ^ 8.0 x 10 [H'^CO+] ^ 6.5 x 10 " and thus 
a [DCO^]/[HCO^]-0.02 abundance ratio towards the core peak. 
Such gas-phase DCO^ enhancement is achieved via reaction: 



Ht H- HD ^ H2D+ -I- H2 + AE 



(1) 



which is endothermic by ~232K in the right-to-left direction 
jGerlich et al. 20021 ), followed by 



CO + H2D+ ^ DCO+ + H2 



(2) 



which dominates the DCO^ formation in the cold and dense gas. 
The absence of significant DCO^ line emission in the PDR is 
consistent with the higher temperatures (>60 K) expected in the 
illuminated edge of the cloud. 

Besides, the detection of intense H'^'CO^ emission towards 
the shielded core and its vicinity (see Fig. [U implies low 



- PDR-like environments require to consider inelastic collisions with 
Hj, H, He and e . H'^CO^, DCO+ and HOC+ collisional rates with H2, 
H and He have been scaled from those of Flower (1999), while colli- 
sional rates with e- were kindly provided by A. Faure and J. Tennyson 
(see e.g., [Faure & Tennyson 200T) . 
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ionization fractions. In terms of excitation and opacity ef- 
fects, H'-^CO^ is a much more rehable tracer of HCO^ col- 
umn density than H'^CO^ itself (as the latter suffers from very 
large opacities and line photon scattering by low-density halos; 
e.g., ICernicharo & Guehn 1987| l. In terms of its chemistry, two 
main processes dominate the formation of H'^CO^ in the low 
temperature shielded gas: 



'COh-H 



+ H2 



and isotopic fractionation through 



(3) 



(4) 



which is endothermic by only ~9K in the right-to-left direc- 
tion (Langer et al. 1984) and competes with dissociative recom- 
bination in the destruction of H'^CO^ where the abundance of 
electrons is low. For the physical conditions prevailing in the 
shielded core, we predict [H'^C0^]/[H'^^C0^] abundance ratios 
down to ~0.7 times lower than the elemental ['^C]/['^C] iso- 
topic ratio. Since both H'^CO^ and DCO^ are mainly destroyed 
by fast dissociative recombination with electrons: 



H'^CO+ + e 



13 



CO-H H 



DCO+ + e" ^ CO + D 



(5) 
(6) 



their abundances inversely scale with that of electrons. In this 
work we have used a "standard" HCO^ dissociative recom- 
bination rate (i.e., a(HCO+)=2.4xlO-^(300/T)°<''* cm^ s"') rec- 
ommended in most astrochemical data bases. We note, how- 
ever, that there is a certain discrepancy among different theo- 
retical calculations and measurements of this key chemical rate 
(see discussion by Florescu-Mitchell & Mitchell 2006 and refer- 
ences therein). In Sect. [T] we discuss the influence of adopting a 
smaller, "non standard" a'(HCO^) rate on our results. 

4.2. HOC+ and H^^CO^ chemistry in fhe PDR 

In order to extract the [HOC+] and [H'^CO+] abundances 
towards the Horsehead PDR, we have modeled the ob- 
served lines (Fig. |2]i using our best knowledge of the pre- 
vailing physical conditions: Tt=60-120K, n(H2)=5xlO'*cm"^, 
n(H)=500cm"^, [e"]=5xlO"^ and a 0.1 pc line-of-sight depth 
(or A^// ^3.1x10^^ cm"^) all accurate within a factor ~2. From 
the observed lines we derive the following column den- 
sities: A^(HOC+)=(1.2-2.5)xlO"cm-2 and A?(Hi^CO+)=(4.7- 
7.8)xl0" cm-2, which translates into [HOC+]=(0.4-0.8)xlO-" 
and [H'^CO+]=(1.5-2.5)xlO"". This computation assumes that 
the HOC+ and H'^CO+ emission fifls the IRAM-30m beam. 
However, HOC^ has not been mapped and its emission could 
well arise from the same ~12"-width filament where the 
emission of small hydrocarbons and HCO radical is concen- 
trated ( jPety et al. 2005| iGerin'et al. 2009l l. In this case, [HOC+] 
increases by a factor ~3. Therefore, we conclude that the 
[HOC-^l/IH'^CO+l abundance ratio towards the PDR lies in the 
range ^0.3-0.8. These values are orders of magnitude higher 
than the value expected in the UV shielded gas. 

Our chemical models (see next section) reproduce the 
[HOC+]/[H'^CO+] abundance ratio towards the PDR but the 
absolute abundances derived from observations are larger than 
those predicted by the model. The discrepancies between ob- 
served and modeled abundances for HOC^ and H'^CO^ have 
likely a common origin. In particular, the formation of HOC^ in 
UV irradiated gas is driven by reactions involving C^ and species 



such as H2O and CO^ (from C^ + OH reaction) that efficiently 
form at high temperatures, that is: 

C+ + H2O ^ HCO+/HOC+ + H (7) 

C0+ + H2 ^ HCO+/HOC+ + H (8) 

where reaction [T] predominantly produces HOC^ whereas reac- 
tion [8] has similar branching ratios for the HCO^ and HOC^ 
formation (e.g., IScottetal. 19971 [Savage & Ziurys 2004| l. The 
HOC^ destruction is governed by the isomerization reaction: 



HOC+ -i- H2 ^ HCO+ + H2 



(9) 



Laboratory experiments show that the reaction rate is smaller 
than previously thought (ISmith et al. 2002l l. allowing interstellar 
HOC^ to exist at detectable amounts. 

The intensity peak of the CO 7=2-1 optically thick lines ob- 
served with the PdBI (T,„i^60K<T„;|Petyetar2005'), together 
with the observed CO 7=4-3/2-1 line ratio (Philipp et al. 2006), 
provide a lower limit to the gas temperature in the PDR (T^ ^60- 
120 K). Temperatures in this range are predicted by the PDR 
model but are not enough to overcome the activation energy 
barriers of the neutral-neutral reactions leading to the forma- 
tion of abundant H2O OH and CO+ (e.g., INeufeld et al. 19951 
ICernicharo et al. 2006I ). Therefore, our models predict HOC^ 
and H'^CO"^ abundances lower than observed because their pre- 
cursor molecules have low abundances, and reactions|2]and[8]are 
not efficient enough. 

We have computed that gas temperatures around ~350 K are 
needed to reproduce the observed HOC^ and H'^CO^ abun- 
dances in the PDR through the previous scheme (see Fig. [3]). Our 
models of the Horsehead (low UV radiation field) include photo- 
electric heating from PAHs and grains but do not predict such a 
warm gas component even if the PAH abundance is significantly 
increased. However, we do not model the PDR gas dynamics and 
thus processes such as shock waves, driven by the expansion of 
the H II region that compress the cloud edge, may provide addi- 
tional gas heating sources to trigger this warm chemistry. This 
reasoning is partially consistent with the non detection of CO^ 
lines, at least at the sensitivity level of our long integration ob- 
servation (rms~50mK in a 0.20 kms"' velocity width channel 
or[CO+]<5xlO-i^). 



H CO and HOC in the PDR 
h"cO* (x60). 



- observed 



X 10 
a 




gas temperature (K) 

Fig. 3. Predicted H'^CO+, HOC+ and CO+ peak abundances in 
the PDR (Av ^0.5-1 .5) as a function of T^,„. H'^CO^ and HOC+ 
abundances (and CO^ abundance upper limit) derived from ob- 
servations towards the PDR are shown with horizontal thin lines. 



If the gas in the PDR has not gone through such a warm 
phase, reaction [8] has to be ruled-out as the main chemical path- 
way for HOC^ formation and an alternative formation scenario 
is required. In this case, we propose that the enhanced HOC^ 
abundance in the PDR can still be related with the high abun- 
dance of C^ (and thus high ionization fraction), but also with 
grain photodesorption of water-ice mantles formed in earlier 
evolutionary stages of the cloud. 
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Fig. 4. Grid of chemical models for 2 different metal abundances, low-metallicity with [M]=10-^ (strong metal depletion case, 
solid curves) and high-metallicity with [M]=10-^ {weak metal depletion case, dashed curves). "Low CR rate" refers to models with 
^=3x10 '^ s ' while models with "high CR rate" refer to ^=3x10 '^ s '. Models "with PAH" include PAH-gas interactions in the 
chemical network (with [PAH]=10-^) while models with "no PAH" are pure gas-phase models. The black empty square represents 
the [H0C^]/[H'''C0^] abundance ratio infeiTed towards the PDR from observations. The blue and red empty squares represent the 
DCO"^ and H'^'CO^ abundances derived towards the shielded core. 
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In this picture, the low xl^ ratio in the Horsehead 
(~10"^) will allow water-ices to be photodesorbed close 
to the illuminated edge of the cloud (see predictions by 
iHollenb ach et al. 2 009 ), increasing the water vapor abundance 
well above the pure-gas phase predictions. Reaction |7] will then 
dominate the HOC^ formation in the PDR. Taking into account 
that isomerization, dissociative recombination and photodisso- 
ciation contributes to HOC^ destruction, we estimate that the 
required water vapor abundance needed to explain the inferred 
HOC+ abundance in the PDR is [H2O]^1800x[HOC+]^(0.7- 
2.2)xl0"^. Herschel Space Observatory observations will 
enable the detection of and H2O lines in a large sample of 
PDRs, confirming whether or not water vapor is abundant at the 
edges of molecular clouds (e.g.. lCernicharo et al. 2006] l. 



5. Determination of the ionization fraction 

Figure |4]presents depth-dependent predictions of several photo- 
chemical models across the Horsehead edge. Each model shows 
the main physical parameters (density and temperature), the ion- 
ization fraction gradient, the DCO+, H'^CO+ and HOC+ abun- 
dances (our observational probes of the ionization fraction) and 
the abundances of key chemical species for the charge balance 
in the cloud: CO/C/C+, M/M+, PAH7PAH/PAH+. 

Four sets of models are displayed. Top/bottom models use a 
low (f =3x10"'^ s"') and high (^=3x10"'^ s"') cosmic-rays ion- 
ization rate respectively. Left/right models exclude and include 
the effects of PAHs respectively. In the latter case, we include 
PAHs in the UV radiative transfer (as a source of absorption 
and scattering of UV photons), in the photoelectric heating and 
in the chemical network. We start the chemistry computation 
by including neutral PAH alone with an initial abundance of 
[PAH] = 1 0"^ . In each set of models (each panel), the only param- 
eter that varies is the abundance of metals: high metallicity with 
[M] = 10"'' (dashed curves) and low metalhcity with [M] = 10"^ 
(solid curves). The low metallicity case implies a large metal 
depletion from the gas phase. 

In terms of the chemical species observed in this work, 
a salient feature of all models is the constancy of the 
DCO^/H'^CO^ abundance ratio once the gas is shielded from 
UV radiation (Ay >6). This feature agrees with the almost identi- 
cal spatial distribution of DCO^ and H'^^CO^ emission observed 
beyond the PDR (see Fig.!!). This similarity was already noticed 
in the lower resolution DCO^ and H'^CO^ pioneering maps of 
several molecular clouds (e.g., IGuelin et al. 1982b . Besides, the 
predicted [HOC+]/[H^^CO+] abundance ratio towards the PDR 
is in good agreement with the value inferred from observations. 
In this UV irradiated region where the and electron abun- 
dances are very high, the HCO^ destruction rate becomes com- 
parable to the isomerization rate (reaction |9]l. This implies that 
the [H0C^]/[H'^C0^] abundance ratio in the cloud achieves the 
highest value in the PDR. 

5.1. The role of ionized carbon and metals 

According to the ionization fraction gradient all models show 
two differentiated environments separated by a transition region: 
the "PDR" (Av -0-2) and the "shielded core" (Ay >6). The elec- 
tron density at every cloud position is given by the difference of 
cations and anions densities; 



In the PDR, carbon, the most abundant heavy element with a 
ionization potential below 13.6 eV, provides most of the charge: 
n(e ) ^ «(C^). Therefore, the ionization fraction in the PDR is 
high, [e ]~10 (or n(,~l-5 cm"^^), and independent of elemental 
abundances other than carbon. 

As Av increases inwards the cloud, the abundance de- 
creases by several orders of magnitude and so does the abun- 
dance of electrons. In the shielded core (Av >6), low ioniza- 
tion heavy metal ions (e.g., Fe"*", Mg'*" or Na^) determine much 
of the ionization fraction (Oppenheimer & Dalgarno 1974 
IGuelin et al. 19821 . In the absence of PAHs, abundant molecu- 
lar ions transfer charge fast to heavy metal atoms M through 

+ M ^ m + reactions. Metal ions recombine orders of 
magnitude slower than molecular ions (TableUJi, and thus a large 
fraction of them is kept ionized (higher [M] implies higher elec- 
tron abundances). Therefore, the ionization fraction in the core 
is highly dependent on the adopted metallicity, and varies from 
a few xlO-'' for [MJ^IO"'', to a few xlO"^ for [MJ^IO^^ 

5.2. The role of PAHs 

Depending on their abundances, the presence of PAHs can al- 
ter the chemistry and the ionization balance in dense clouds 
(e.g., |Lepp & Dalgarno 1988| l. For our adopted abundance of 
[PAH] = 10 ' the right and left panels in Fig. |4] shows that the 
presence of PAHs mostly modifies the ionization fraction at 
Av >2. If not all PAHs accrete onto bigger grains or coagulate to- 
wards cloud interiors, PAH" can be abundant through the cloud 
because the radiative electron attachment rate 



PAH + e- 



PAH + hv 



(11) 



^ n,(cations^) - ^ nj(anions ). 



(10) 



is high (>10 ^cm^s '), although probably dependent 
on the PAH si ze (Omont 1986' 'Allama ndola et al. 19891 
Flo wer et al. 20071 |mkelam & Herbst 2008). In the shielded 
core PAH" is destroyed by recombination with atomic (M^,...) 
and molecular cations (HCO^, HsO^,...) which are orders of 
magnitude less abundant than the available cations in the PDR 
(C^, S^,...). Negative PAH ions thus reach high abundances 
([PAH"]^2xlO"'^). For our choice of PAH parameters, this 
means that PAH can be the most abundant negatively 
charged species, more than electrons for Av >5. In addition, 
recombination of atomic ions on PAH" is by far more efficient 
than the slow radiative recombination on electrons. This is a 
very important point since heavy metal ions are now neutral- 
ized at similar rates than molecular ions. As a result, both the 
abundance of metal cations and the ionization fraction decreases 
when PAHs are included, while molecular ions such as H'^CO^ 
and DCO^ increase their abundances (see Fig.^right panels). 

In the illuminated edge of the cloud, PAH" is predominantly 
destroyed by UV photons through electron detachment, 

PAH" +hv^ PAH -He" (12) 

and through recombination with atomic cations which are 
very abundant in the PDR (e.g., Bakes & Tielens 19981 
IWolfire et al. 20081 ). In consequence, the abundance of ions such 
as in the PDR decreases with respect to models without 
PAHs. This effect is important to determine the elemental abun- 
dances and depletion factors. Despite the higher PAH" destruc- 
tion rates in the PDR, the high electron density and relatively 
low UV field in the Horsehead allows PAH" to form efficiently 
through electron attachment (reaction fTTTi. Hence the resulting 
PAH" abundance is also high in the PDR. On the other hand, the 
predicted abundance of positively charged PAH^ in our mod- 
els is ~500 times smaller than the abundance of PAH" (due to 
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Fig. 5. Left: H CO^ and DCO^ spectra along the direction of the exciting star at 6y=l5" (histograms). Radiative transfer models 
using the output of PDR models for a fixed cosmic rays ionization rate (^=3x10"'^ s"') and varying metallicities. Thin blue curves 
for [M]=5xlO"^ and no PAHs; thick grey curves for [M]=10"'' and no PAHs; and dashed red curves for [M] = 10"'' and [PAH] = 10"^. 
Modeled line profiles have been convolved with the appropriate Gaussian beam at each observed frequency (the angular resolution 
for each line are quoted in Tables[T]and|2]l. Right: Same as previous figure but for a fixed metal abundance ([M]=lxlO"^), no PAHs 
and varying cosmic rays ionization rate Thin blue curves for a model with ^=5x10""^ s"'; thick grey curves for ^=3x10"'^ s"'; 
and dashed red curves for ^=10"'^ s"'. 



fast electronic recombination) and therefore PAH^ don't seem to 
play a major role in the ionization balance inside the cloud (see 
also|Lepp & Dalgarno 198"8llWakelam & Herbst 2008]l. 



5.3. The role of the cosmic-ray ionization rate 

Cosmic rays affect the ionization state and the physics of molec- 
ular clouds, being the dominant source of heating and ioniza- 
tion in the gas shielded from interstellar radiation fields. Indeed, 
secondary UV photons are created in cloud interiors by H2 
electronic cascades following H2 excitation by collisions with 
cosmic rays (IPrasad & Tarafdar 1983 1. Therefore, cosmic rays 
maintain a certain ionization degree in the shielded gas and play 
a fundamental role in the ion-neutral chemistry by setting the 
abundance of key ions ( [Herbst & Klemperer 1973 1. 

Most studies based on the interpretation of observed 
molecular ions set a range of a few 10"^'^ to a few 
10 s ' for the cosmic -ray ionization rate (|Le Petit et al. 20041 



Ivan der Tak 20061 [Dalgarno 2006 and references therein) 



However, it is still discussed whether or not ( depends on envi- 
ronmental conditions {e.g., galactic center vs. disk clouds) or if 
it varies from source to source (e.g., from dense molecular cores 
to more translucent clouds). In many ways, PDRs offer an in- 
teresting intermediate medium to analyze the transition between 
translucent and dark clouds. 

In terms of our observations, the DCO^ and H'^CO^ abun- 
dances directly scale with ^ in the UV shielded gas. Indeed, these 
ions are direct products of the Hj destruction (through reac- 
tions |2] and O, and the HJ formation is proportional to mh- 
However, ^ and the metal abundance can not be constrained in- 
dependently from the infeiTed DCO^ and H'^CO^ abundances 
since both parameters control the ionization fraction, and thus 
the destruction of these ions through reactions|5]and|6] 

6. Results: observational constraints 

In this section we compare the synthetic and observed H'-'CO^ 
and DCO^ spectra as a function of cloud position. We then ex- 



plore the range of metallicities and cosmic -rays ionization rates 
compatible with the H'^^CO^ and DCO^ inferred abundances 
(see Table|6l). The influence of PAHs is also investigated. We fi- 
nally compare the [HOC^]/[H'^CO^] ratio obtained towards the 
Horsehead with the values derived in other PDRs. 



6. 1 . Constraints to the metals abundance 

Figure |5] left shows the spectra along the direction of the ex- 
citing star (histograms) and radiative transfer models using 
the output of several PDR models for a fixed ionization rate 
(^=3x10"'^ s"') and varying metallicities. In particular, the 
model with [M] = 10"' (and no PAHs) displays a notable agree- 
ment with both the DCO^ and H'^CO^ spatial distribution and 
with the inferred peak abundances towards the core (Table |6]l. 
In addition. Fig. |6] left shows the predicted ionization frac- 
tion and [H'^CO^] and [DCO^] abundances at the core peak 
(Ay >10) as a function of [M] (blue-solid curves). These mod- 
els (no PAHs, fixed f) shows that the upper limit metallicity 
compatible with observations is [M]< 4x10^, which implies 
a ionization fraction of [e ] =(7+l)xl0^ at the core peak. 
Higher metal abundances increase the ionization fraction (see 
Fig. |6] left), which translates into weaker lines than observed 
(Fig.|5]fe/f: blue-thin curves). Therefore, the gas-phase metal 
abundance is depleted by ~4 orders of magnitude with respect to 
the Sun ([M] ^ 8.5x10'^; [Anders & Grevesse 1989l l. This range 
of depletion is similar to that obtained in other prestellar cores 



such as Barnard 68 (Maret & Bergin 2007 1. We shall refer it as 
the strong metal depletion case. 

The inclusion of PAH interactions implies lower ioniza- 
tion fractions and enhanced molecular ion abundances (see 
Fig. [6]Ze/0 which results in overestimated H'^CO^ and DCO^ 
line intensities towards the core. Therefore, the abundance of 
metals {e.g., the ionization fraction indirectly) has to be in- 
creased to match the observed intensities. In particular, Fig.[5]fe/f 
shows that a model with [PAHs] = 10"^ and [M] = 10"'' (red- 
dashed curves) displays only a factor <2 brighter lines than 
models with 3 orders of magnitude lower metallicities and no 
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Fig. 6. Left: Model predictions for the shielded core (the "DCO^ peak" at Ay > 10). The ionization rate due to cosmic rays is fixed to 
^=3x10 s ' . The different panels show {upper): the ionization fraction, (middle): the H'^CO^ abundance and {lower): the DCO^ 
abundance as a function of gas phase metallicity. Blue-solid curves for models without PAHs, and red-dotted curves for models 
with neutral and charged PAHs ([PAHJ^IO"^). Horizontal shaded regions show the H'^CO^ and DCO^ abundances derived from 
observations towards the core, while vertical shaded regions show the parameter space compatible with observations. Right: Same 
as previous figure but for a fixed low metallicity of [MJ^IO"'' (no PAH; blue-solid curves) and a fixed high metallicity of [M] = 10-* 
([PAH] = 10"^; red-dotted curves). The different panels show {upper): the ionization fraction, {middle): the H'^CO^ abundance and 
{lower): the DCO^ abundance as a function of the ionization rate due to cosmic rays. 



Table 6. Infen-ed 
Nn^N{Yi)+2N{}\2). 



abundances [x] =A^(x)/A^h where 



Species 


shielded core 


PDR 




Av >6 


Av=0-2 




Sx^AS" 


6x^\5" 


Nh (cm--) 


5.8 X 10^- 


3.1 X 10^- 


[H"CO+] 


6.5 X 10-" 


1.5 X 10-" 


[H'^CO+l 


3.9 X 10-' 


9.0 X 10-'" 


[DCO+] 


8.0 X 10-" 


(-) 


[HOC+] 


(-) 


0.4 X 10-" t 


[CO+] 


(-) 


< 5.0 X 10-" 




(1 - 8) X 10"'' 


10-" - 10-^ 



t Assuming extended emission. It would be 1.2x10 " if HOC* 
from a 12"-width filament as HCO dGerin et al. 2009l l. 



PAHs (grey-thick curves). As shown in Fig. |6] left, the inclu- 
sion of PAHs makes the range of metal abundances compat- 
ible with observations much higher now, [M]=(3+l)xl0^. 

The required abundance of metals is at least a factor ~500 larger 
than the former case without PAHs. Thus, we refer it as the weak 
metal depletion case. Note that this metallicity is still below the 
gas-phase abundance of Fe-i-MgH-... elements in the diffuse in- 
terstellar gas {e.g., IWolfire etal. 19951 IHowk et al. 20061 1 and is 
consistent with the incorporation of metals into dust grains in 
higher density regions {e.g., IWolfire et al. 199 5). On the other 
hand, the ionization fraction required to reproduce the H'^CO* 
andDCO* abundances does not change much, [e-]=(4±l)xl0-'' 
at the core peak (red-dashed curves Fig.|6]Ze/0- Therefore, it is 
not easy to distinguish between the strong metal depletion (no 
PAHs) and weak metal depletion (with PAHs) cases in terms 
of the ionization fraction. The observation of forbidden lines 
from metals such as [Fe ii] towards the UV illuminated edges of 



molecular clouds may help to remove this apparent degeneracy. 
In one of the few positive cases, the S140 PDR, the detection 
of a weak [Feii]26.0yum fine-structure line emission suggests 
that iron is depleted, but with an abundance of -5x10-*^ relative 
to H (ITimmermann et al. 19961 ). Nevertheless, without mapping 
and comparing with other PDR tracers, it is not obvious to dis- 
entangle whether [Fe ii] lines arise from the PDR gas or from 
the adjacent (Hn) ionization front {e.g., IMarconi et al. 199"8l 
IKaufman et al. 2006l l. 



6.2. Constraints to the cosmic-ray ionization rate 

Figure^ right shows again the observed spectra along the di- 
rection of the exciting star (histograms) and radiative trans- 
fer models using the output of PDR models that vary the cos- 
mic rays rate, without PAHs and a fixed metal abundance of 
[M] = 10-^. The adopted metal abundance is, within an order 
of magnitude, the usual value estimated in prestellar cores 
{e.g., ICaselli et al. 19991 |Maret & Bergin 2007[ ) and is compat- 
ible with our strong metal depletion case. Therefore, Fig.^right 
shows the effects of different ionization rates directly on the 
DCO* and H'^CO* line intensities. For the adopted physical 
conditions and chemical network, our observations constrain ^ 
within a factor ~2. In particular. Fig. |6] right shows the pre- 
dicted [e-], [H'^CO*] and [DCO*] abundances at the core peak 
(Ay > 10) as a function of and evidences that in the absence of 
PAHs, the cosmic rays ionization rate range compatible with 
the observations of the Horsehead edge is ^=(5+3)xl0 s ' 
(blue-solid curves). If PAHs are included in the chemistry, the 
metal abundance has to be increased accordingly to reproduce 
the observations. For our [PAH] = 10-^ model case, the required 
metal abundance needs to be above [M]^10-^ to obtain ioniza- 
tion rates above ^^10-' 



{Fig-Aright; red-dotted curves). 
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Note that given the fact that the H'^CO^ formation in the 
PDR is not dominated by the '^CO + H+ reaction, the H'^CO+ 
abundance in the UV illuminated gas does not scale with ^. 
Therefore, we can not further investigate if f varies signifi- 
cantly in the transition from diffuse regions to the shielded core 
(e.g- IMcCall et al. 2003llPadoan & ScalolOOSI) . 

6.3. High ionization fraction in fhe PDR 



The bright [Cii]15 8/jm dZhou et al. 199? ) and [Ci]492GHz 
( [Philipp et al. 2006| fine structure line emission towards the 
Horsehead PDR, together with subtle chemical effects such 
as the large [/-C3H2]/[c-C3H2] linear-to-cyclic abundance ratio 
( [Teyssier et al. 2005) all show observationally that the ionization 
fraction is higher in the UV illuminated edge than towards the 
cloud interior. Nevertheless, all those studies lacked the angular 
resolution to properly measure the ionization fraction gradient. 

The abundances of reactive ions such as HOC^ are also 
predicted to be enhanced in the UV illuminated gas, where 
we have shown that the ionization fraction is high, up to 
[e ]~10"^, and that the HOC^ formation is linked to the avail- 
ability of C^. On the other hand, the H'^CO^ abundance in- 
creases as the electron abundance decreases towards the shielded 
core. Therefore, we predict that the [HOC+]/[H'^CO+] abun- 
dance ratio scales with the ionization fraction gradient, 
reaching the highest values in the PDR (Fig. |4|. In particular, 
we derive a high [HOC+]/[H^^CO+]=0.3-0.8 ratio (or a low 
[HCO+]/[HOC+] -75-200 ratio) towards the PDR, similar to tha t 
observed in other PDRs such as NGC7023 jFuente et al. 20031 1. 

7. Discussion 

7.1. Tlie ionization fraction gradient 

Star forming clouds display different environments depending 
on the dominant physical and chemical processes. These envi- 
ronments are, in a first approximation, similar to those studied 
here: (/) a low density cloud edge directly exposed to a UV ra- 
diation field from nearby stars; (//) a transition region or ridge 
where the H2 density increases as the gas temperature decreases 
due to the attenuation of the external radiation field. UV pho- 
tons can still play a significant role depending on their penetra- 
tion depths (e.g., cloud dumpiness, grain properties, etc.) and 
(///) the denser shielded cores that may be externally triggered 
to form a new generation of stars depending on their stability 
against gravitational collapse (e.g.. |Goicoechea et al. 2008] l. 

Assuming that the observed field-of-view in the Horsehead 
nebula is representative of the above 3 environments, our maps 
and chemical models reveal that the ionization fraction fol- 
lows a steep gradient in molecular clouds: from [e ] ^10 
at the edge of the cloud (the "C^ dominated" region) to a few 
times ~10"^ in the shielded cores. The prevailing chemistry 
and the abundance of atomic ions such as and deter- 
mine the "slope" of the ionization fraction gradient in the tran- 
sition regions. In particular, sulfur (with a ionization poten- 
tial of ~ 10.36 eV) is a good source of charge behind the "C^ 
dominated" region. Advection and time-dependent effects may 
also modify the ionization fraction gradient with time. However, 
Morata & Herbst (2008) have shown models for (uniform) phys- 
ical conditions similar to those in the Horsehead where [HCO^] 
(and [e^] presumably) does not change with time appreciably. 

Figure [T] shows the predicted ionization fraction gradi- 
ent, with a scale length of ~0.05pc (or ~25"), and the main 
charge carriers for 3 representative models with fixed stan- 



dard metal abundance (our strong metal depletion case) and 
standard cosmic -ray rate. Panel [TJj shows a model without 
PAHs and high gas-phase sulfur abundance ([S]=3.5xlO"''; 
IGoicoechea et al. 20061 1. The ionization fraction gradient in 
the core / transition / edge zones, is mainly determined by 
the [HCO+-hS+-hM+-h...]/[S+]/[C+] abundances respectively. 
Sulfur ions control the charge balance in the transition layers, 
and due to their high abundance and slow radiative recombina- 
tion rate with electrons, the ionization fraction is high, a few 
times 10"^, and the gradient is smooth. This model qualitatively 
agrees with the observed more extended emission and narrower 
line-widths of sulfur recombination lines compared to carbon 
recombination lines in dark clouds, as well as with the rela- 
tively low (< 10) carbon-to-sulfur recombination lines intensity 
ratio (e.g., Pankonin & Walmsley 1978 Falgarone et al. 1978 1. 
All these signatures argue in favor of extended regions with sig- 
nificant amounts of gas-phase S^. 

Panel [Tj? shows the same model but sulfur abundances 
smaller by two orders of magnitude (strong gas-phase sulfur 
depletion). The ionization fraction in the cloud core (Ay >6) 
and edge (Ay <2) are nearly the same as in the previous high 
sulfur abundance model. However, the lack of abundant S^ in 
the transition layers decreases the electron abundance consider- 
ably, and makes the ionization fraction gradient much steeper 
Observational constraints to the atomic and ionic S abundances 
from far-IR fine structure lines or recombination lines, and a 
careful treatment of the sulfur chemistry (i.e., which are the most 
abundant S-bearing species as a function of cloud depth?) are 
thus required to quantify the S^ abundance at large Ay and its 
impact on the charge balance. 

Panel |7]: finally shows a model with high sulfur abun- 
dance again but including PAHs (with [PAH] = 10"^). As pre- 
sented in Sect. 15.21 negatively charged PAH" efficiently form 
by radiative electron attachment and their abundance remains 
high through the cloud. Given the much higher recombination 
rates of atomic ions on PAH" than on electrons, the abun- 
dance of atomic ions such as S^ in the transition zone, or 

in the shielded cores, quickly decreases. Hence, lower ion- 
ization fractions (and a much weaker dependence on the as- 
sumed metal elemental abundance) are predicted by the model 
with PAHs. This results agrees with theoretical predictions for 
UV shielded gas (|Lepp & Dalgarno 19881 IFlower et al. 2007] 
IWakelam & Herbst 2008l l. 

In summary, a high abundance of PAHs throughout the 
molecular cloud (not only in the PDR) plays a role in the ion- 
ization balance and in the abundance of molecular ions, which 
affects the determination of elemental abundances (e.g., S) from 
fractional molecular abundances (e.g., HCS^/CS, SOVSO, etc.). 

7.2. Tile PAH abundance in UV shieided gas 

The PAH abundance in the dense and UV shielded gas is 
far from being well constrained. Different approaches to an- 
alyze ISO and Spitzer mid-IR observations towards several 
PDRs all argue in favor of an evolution of dust grain sizes: 
from the illuminated cloud edge where the emission of PAH 
bands dominates, to the shielded interiors where the continuum 
emission from bigger grains dominates ( Rapacioli et al. 2005| 
IBerne et al. 2007. Compiegne et al. 2008). It is not trivial to 
disentangle whether this is a physical effect (i.e., free PAHs 
are not present in the shielded regions) or an excitation effect 
(i.e., lack of UV photons). Even if the PAH abundance drasti- 
cally decrease towards cloud interiors, a chemically significant 
fraction of them may survive. Unfortunately, while the effects 
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Fig. 7. Derived abundance profiles for the most significant ions studied in this work for different PAH and sulfur elemental abun- 
dances. An enhanced UV radiation field 60 times the mean ISRF illuminates the cloud from the right. The metal abundance 
([M] = 10"^) and the cosmic-ray rate (^=3x10"'^ s"') are fixed in all models, a) Model with [PAH]=0 and [S]=3.5xlO"'' (low sulfur 
depletion), b) Model with [PAH]=0 and [S]=3.5xlO-^ (high sulfur depletion), c) Model with [PAH] = 10-^ and [S]=3.5xlO ^ 



of grain growth in the UV extinction curve have been modelled 
by us (Go icoechea & Le Bourlot 20 07). including PAH coagula- 
tion/accretion in the chemistry is beyond the scope of this work. 
All we can say at this point is that a better description of the 
cloud chemistry may be a decreasing PAH abundance gradient 
or an increasing PAH size distribution towards the cloud interior. 
In any case, we have shown than the presence of free PAHs in 
molecular clouds modifies the prevailing chemistry. As a result, 
the predicted high abundance of PAH" can dominate the recom- 
bination of metal ions and reduce the ionization fraction. 

The presence of abundant free PAHs, or large molecules to 
which electron attach (Lepp & Dalgarno 1988), can thus be cru- 
cial in determining the coupling of the gas with magnetic fields 
in molecular clouds, but also in collapsing cores or in the "dead" 
zones of protoplanetary disks (magnetically inactive regions 
where accretion can not occur if the ionization fraction is very 
low). According to our models, the PAH abundance thresh- 
old required to affect the metal and electron abundance de- 
termination in the UV shielded gas is [PAH]> 10^. Herschel 
observations might allow the identification of specific PAH 
carriers through their far-IR skeletal modes dJoblin et al. 20021 
IMulas et al. 2006l l, thus providing clues on their composition 
and abundance variations in different environments. 



7.3. "Non standard" HCO^ dissociative recombination rate 

We conclude by discussing the sensitivity of our determina- 
tion of the ionization fraction from H'^CO^ and DCO^ abun- 
dances. In particular, we have checked the effects of adopting 
a smaller, "non standard" HCO^ dissociative recombination 
rate, a'(HCO+)^0.7xlO-^(300/T)° -^°cm3 s'' OSheehan 20001 
IFlorescu-Mit chell & Mitche ll 2006^. For models without PAHs, 
the predicted H^^CO^ and DCO^ abundances increase by a fac- 
tor ~3 with respect to models using the "standard" a'(HCO^) 
rate (Table HI, but the metallicity required to fit the observed 
lines has to be increased to [MJ^SxlO^** and the predicted ion- 
ization fraction increases to [e^J^SxlO^*^ in the core. This value 
should be regarded as the upper limit of our determination. On 
the other hand, the influence of q''(HCO^) in models with PAHs 
is less important. It also requires high metallicities to fit the ob- 



served intensities (weak metal depletion case), but the predicted 
[e ] in the shielded core remains low (below ~10 '^). 



8. Summary and conclusions 

We have presented the first detection of HOC^ reactive ion 
towards the Horsehead PDR. Combined with our previous 
IRAM-PdBI H'3CO+ 7=1-0 jGerin et al. 2009]l and lRAM- 30m 
H'^CO+ and DCO+ higher-7 line maps (|Pety et al. 2007|) we 



performed a detailed analysis of their chemistry, excitation and 
radiative transfer to constrain the ionization fraction as a func- 
tion of cloud position. The observed field contains 3 different 
environments: (/) the UV illuminated cloud edge, (ii) a transi- 
tion region or ridge and (Hi) a dense and cold shielded core. We 
have presented a study of the ionization fraction gradient in the 
above environments, which can be considered as templates for 
most molecular clouds. Our main conclusions are the following: 

1 . The ionization fraction follows a steep gradient, with a scale 
length of -0.05 pc (-25"), from [e^] ^lO""* («,~l-5 cm^^) at 
the cloud edge (the "C^ dominated" regions) to a few times 
~10"^ in the shielded core (with ongoing deuterium frac- 
tionation). Sulfur, metal and PAH ions play a key role in the 
charge balance at different cloud depths. 

2. The detection of HOC+ towards the PDR, and the high 
[HOC+]/[H'^CO+]-0.3-0.8 abundance ratio inferred, proves 
the high ionization fraction in the UV irradiated gas. 
However, the H'^CO^ and HOC^ abundances derived from 
observations are larger than the PDR model predictions. We 
propose that either the gas is/was warmer than predicted or 
that significant water ice-mantle photodesorption is taking 
place and HOC^ is mainly formed by the + H2O reaction. 

3. The ionization fraction in the shielded core depends on the 
metal abundance and on the cosmic-rays ionization rate. 
Assuming a standard rate of ^=3x10"'^ s"' and pure gas- 
phase chemistry (no PAHs), the metal abundance has to be 
lower than 4xl0~^ (strong metal depletion). Conversely, as- 
suming a standard metal abundance of [M] = 10"^, our obser- 
vations can only be reproduced with ^=(5+3)xlO"'^ s"^ 

4. The inclusion of PAHs modifies the ionization fraction gradi- 
ent and decreases the metal depletion required to reproduce 
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the observations if [PAH]> 10 (i.e., if not all PAHs coagu- 
late/accrete onto bigger grains). In such a case, PAH" acquire 
large abundances also in the shielded gas. Recombination of 
atomic ions on PAH" is much more efficient than on elec- 
trons and thus metal ions are neutralized at similar rates than 
molecular ions. For [PAH] = 10"^, the metal abundance con- 
sistent with observations increases to [M]=(3±l)xlO"^ (still 
below the heavy metals abundance in the diffuse ISM). 
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